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ABSTRACT 

Context. The study of prestellar cores is essential to understanding the initial stages of star formation. With Herschel more cold 
clumps have been detected than ever before. For this study we have selected 21 cold clumps from 20 Herschel fields observed as a 
follow-up on original Planck detections. We have observed these clumps in 13 CO (1-0), C 18 0 (1-0), and NoH + (1-0) lines. 

Aims. Our aim is to find out if these cold clumps are prestellar. We have examined to what extent independent analysis of the dust and 
the molecular lines lead to similar conclusions about the masses of these objects. 

Methods. We calculate the clump masses and densities from the dust continuum and molecular line observations and compare these 
to each other and to the virial and Bonnor-Ebert masses calculated for each clump. Finally we examine two of the fields with radiative 
transfer models to estimate CO abundances. 

Results. When excitation temperatures could be estimated, the column densities derived from molecular line observations were 
comparable to those from dust continuum data. The median column density estimates are 4.2xl0 21 cnT 2 and 5.5xl0 21 cnT 2 for the 
line and dust emission data, respectively. The calculated abundances, column densities, volume densities, and masses all have large 
uncertainties and one must be careful when drawing conclusions. Abundance of 13 CO was found in modeling the two clumps in the 
field G131.65+9.75 to be close to the usual value of 10~ 6 . The abundance ratio of 13 CO and C 18 0 was ~10. Molecular abundances 
could only be estimated with modeling, relying on dust column density data. 

Conclusions. The results indicate that most cold clumps, even those with dust color temperatures close to 11 K, are not necessarily 
prestellar. 
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1. Introduction 


The details of the initial stages of star formation are still poorly 
known, despite many observational and theoretical studies. This 
is in part because star formation involves a very complex mixture 
of effects of gravity, turbulence, rotation, radiation, thermody¬ 
namics, and magnetic fields. Both large representative samples 
of clumps and detailed studies of their properties and relations 
to their cloud environment are needed. 


Herschel observations of star- 
conducted within the Gould Belt 

forming clouds, e.g., those 

Survey (Andre et al. 

2010) 

and HOBYS program (Motte et al. 

[2010), have lead to t 

te de- 

tections of hundreds of starless anc 
This is particularly true for nearby 

protostellar condensations, 
star forming clouds where 


the Herschel resolution is sufficient to resolve gravitationally 
bound cores. The objects appear preferentially within filaments. 


which are a conspicuous feature of all interstellar clouds (Andre 
et al. |2013j ). The filaments themselves contain both supercritical 
gravitationally bound structures and subcritical filaments that are 
likely to disperse with time ( Arzoumanian et a l.|2013j l. 

The Planck satellite has performed an all-sky survey to 
map anisotropies of the cosmic microwave background ( |Tauber| 
et ah] 20101. At the same time Planck is providing maps of 


thermal dust emission from molecular clouds within the Milky 
Way. From these data more than 10,000 compact cold sources 


* Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 


have been detected and a large fraction of these are believed 
to be prestellar cores or larger structures harboring prestellar 
cores ( Planck Collaboration et al.|201 lb) . The Planck and IRAS 
(100 pm) surveys have been used to compose the Cold Clump 
Catalogue of Planck Objects (C3PO), from which 915 most 
reliable detections were published as the Early Cold Clumps 
Catalogue (353-857 GHz 350 -850 pm), ECC (see, e.g., |Planck 


[Collaboration et al.|201 1c|d ej). The sources were detected with 
the method described in Montier et al.| ( |20 1 0| > and the ini¬ 
tial results on the Planck detections were described in IPlanckl 
[Collaboration et al7| ( |201 le|b] l. The final version of C3PO is cur¬ 
rently in preparation and will provide a global view of the cold 
clump population over the whole sky. 


A number of Planck detections of cold clumps have been 
confirmed with the Herschel Space Observatory observations 
(100-500 pm) to be cold (Td ust ~14 K or below) and also dense. 
For the present study, we selected clumps from fields covered 
in the Herschel open time key program Galactic Cold Cores 
(Juvela et al. |2010| ). The aim of this Herschel program is to ex¬ 
amine a representative cross section of the source population of 
cold clumps observed with Planck and to determine the physi¬ 
cal properties of these clumps. The Herschel results suggest that 


many of the sources are already past the prestellar phase (Juvela 
|et al.||20l0| |201 lj l. In this paper, we refer to sources that are 
gravitationally bound as prestellar objects. If the estimated mass 
exceeds the virial mass, the object is expected to be gravitation¬ 
ally bound. Bonnor-Ebert (BE) spheres are used as an alternative 
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model to recognize prestellar cores (|Andre et al.|2010; Konyves 
et al.|2010) . 

In an earlier study, a sample of 71 fields at distances rang¬ 
ing from ~100 pc to several kiloparsecs were examined. The 
Herschel observations, together with AKARI and WISE in¬ 
frared data, were used to confirm the presence of cold dust. In 
about half of the observed fields, point sources were found in 
the mid-infrared, indicating active star formation. However, the 
cold dust still dominated the submillimeter spectra in these ac¬ 
tive star formation areas ( [Juvela et al.|2012b| . 

Some of the Planck- detected clumps have already been stud¬ 
ied in molecular lines. Wu et al. ( 201 2[ ) surveyed 673 sources 
with single-point observations and mapped 10 clumps with 22 
identified potential cold cores. They found seven cores that are 
likely to be gravitationally bound and, thus, on the verge of col¬ 
lapse. In a follow-up study, 71 of the sources observed only with 
single-point observations were examined and 90 % of the found 
cores were starless (Meng et al.||2013l. The clumps studied by 


Wu et al. (2012l and Meng et al. ( 2013| l did not include the 


clumps chosen for this paper. 

The definition of dense structures in the ISM is still an on¬ 
going process. The clumps are small and dense condensations 
in the molecular clouds, often inside filaments. Cores are even 
smaller and denser. |Bergin & Tafalla ( 2007) 1 defined the clumps 
to have a mass of 50-500 Mq, size of 0.3-3 pc, and density of 
10 3 —10 4 cm~ 3 . Cores were defined to have mass of 0.5-5 M Q , 
size of 0.03-0.2 pc, and density of 10 4 —10 5 cm~ 3 . We use the 
word clump as a general term for the dense structures we inves¬ 
tigate here, however, some of these structures are likely to be 
cores themselves or contain unresolved cores. 

Extensive surveys have been done in 12 CO and l3 CO over 
the Galactic plane and of nearby star-forming clouds (see, e.g., 
Combes[l991 for review). However, general surveys are not the 


best method to research possible compact sources. For the first 
time, with Planck and Herschel , we can identify the coldest 
and densest clumps in the known molecular clouds complexes 
and focus on dedicated molecular line observations targeting the 
most interesting dense clumps in these regions. 

We combine Herschel data with molecular line observations 
to examine the physical properties of several cold clumps and 
their likelihood to evolve into star-forming cores. To investigate 
the gravitational stability of the objects, molecular line obser¬ 
vations are essential to providing direct information of the ki¬ 
netic and turbulent forces and, thus, on the current stability of 
the clumps and the cores. The objects can become gravitation¬ 
ally unstable by cooling and accretion, possibly aided or hin¬ 
dered by turbulence and external forces ( Bergin & Tafalla|2007 
|Ballesteros-Paredes et al.|2007) l. 

We investigate the physical state of a few clumps that were 
originally selected based on the cold dust signature detected by 
the Planck satellite. We derive the CO column densities, vol¬ 
ume densities, and clump masses and estimate the gravitational 
stability. The potential prestellar nature of the objects is inves¬ 
tigated by comparing their mass estimates with virial masses. 
We look for signs of CO abundance variations as further evi¬ 
dence of the nature of the sources. Analysis assuming a local 
thermodynamic equilibrium (LTE) will be complemented with 
modeling. We selected one target for radiative transfer modeling 
where continuum data are used to constrain the models so that 
direct estimates of CO abundance are possible. The results also 
provide insights into the nature of the larger Planck cold clump 
population. 

The paper is organized as follows. The observations are de¬ 
scribed in Sect. [2] Further, we go through the methods in Sect. 


[3] The results are presented in Sect. [4] and they are discussed in 
more detail in Sect. [5] Finally, the conclusions are summarized 
in Sect. [6] 


2. Observations 

2.1. The clump sample 

We examine 21 clumps in 20 fields previously observed with 
Herschel at wavelengths 100 - 500 pm. The fields mapped 
with Herschel were originally selected based on the Planck 
survey, in which these clumps showed a significant excess of 
cold dust emission ( [Planck Collaboration et aT~| [2011 c[ >. Our 
selection was based on the Herschel data from the Galactic 
Cold Cores Herschel key program (PI Juvela) that carried out 
follow-up observations of 116 fields, each with one or more 
Planck cold clumps. We chose the bright clumps, which based 
on Herschel data also contain very cold dust (color temperature 
at Tdust ~14 K, 40" resolution). The fields are located in the 
Milky Way at galactic latitudes |b|=9 - 17°, which ensures min¬ 
imal confusion from background emission. While they typically 
have a general cometary or filamentary morphology, the submil¬ 
limeter data show that at small scales the clouds are fragmented 
to several clumps. 


2.2. Molecular line observations 


We mapped the selected clumps in the 13 CO (1-0) spectral line 
using the 20-m radio telescope in Onsala Space Observatory, 
Sweden, in February and March 2012. Toward the 13 CO peaks, 
we also observed C ls O (1-0) and N 2 H + (1-0) lines. Part of the 
northern clump in field G131.65+9.75 was mapped in 13 CO (1- 
0) in February 2011 (Planck Collaboration et al. 201 lb i. The 
observed clumps are summarized in Table |T] and the observed 
spectra are shown in the appendix [A] We carried out the obser¬ 
vations in frequency switching mode, with frequency throws of 
+ 10 MHz for l3 CO and C lg O and +5 MHz for N 2 H + . For N 2 H + 
in the 86 GHz band the half-power beam width, HPBW, is 44" 
and the main beam effciency, ?/ B , is 0.65. For 13 CO and C 18 0 in 
the 109 GHz band the HPBW is 35" and p B is 0.45. 

The telescope has a pointing accuracy of 3" rms. The 
pointing was checked with SiO spectra toward bright sources 
R Cassiopeiae and U Orionis. The calibration was achieved 
through the chopper-wheel method, and the system gives an¬ 
tenna temperatures in units of T* A (Kutner & Ulich|| 1981 i. The 
backend was a low-resolution digital autocorrelator spectrome¬ 
ter (ACS) with a bandwidth of 40 MHz divided into 1198 chan¬ 
nels. This gave a channel width of 25 kHz, corresponding to 
0.068 km s _I at the frequency of the 13 CO(l - 0) line. 

On average the RMS noise for the observed 13 CO lines was 
~0.27 K per channel on the main beam temperature 0/ \ib) scale. 
The corresponding values for the observed C ls O and N 2 H + 
lines are ~0.06 K and ~0.05 K, respectively. The data were re¬ 
duced with GILD AS/CLASS software package^] The frequency 
switched spectra were folded, averaged, and the spectral base¬ 
lines were modeled and subtracted using third order polynomi¬ 
als, except for C ls O in sources G108.28+16.68, G131.65+9.75 
north and south, and G161.55-9.30 where fifth order polynomi¬ 
als were used. 


1 Institut de Radioastronomie Millimetrique (IRAM): 
http://iram.fr/IRAMFR/GILDAS/ 
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Table 1 . Summary of observations. The table lists the 21 targets, kinematic distances for the fields (Montillaud et al. 2015), temper¬ 
ature of the dust, coordinates for the 13 CO (1-0) peak temperature (where the C ls O (1-0) and N 2 IF (1-0) were observed), area in 
square arc minutes selected for the 13 CO (1-0) mapping, and number of observed points for C 18 0 (1-0) and N 2 H + (1-0). 


Field 

d (pc) 

Tdust 

a (J2000) 

5 (J2000) 

13 CO (sq') 

C 18 0 (N) 

N 2 H + (N) 

G86.97-4.06 

700 

11.5+0.9 

21 17 44.22 

+43 18 47.1 

4.55 

1 

1 

G92.04+3.93 

800 

10.4+0.8 

21 02 23.43 

+52 28 38.3 

3.96 

1 

1 

G93.21+9.55 

440 

10.6+0.8 

20 37 19.98 

+56 54 42.6 

6.79 

2 

1 

G94.15+6.50 

800 

12+1 

20 59 04.80 

+55 34 33.4 

6.67 

1 

1 

G98.00+8.75 

1100 

12+1 

21 04 22.70 

+60 08 55.1 

3.79 

1 

1 

G105.57+10.39 

900 

10.9+0.8 

21 41 01.79 

+66 34 40.9 

8.33 

1 

1 

G108.28+16.68 

300 

14+1 

21 09 08.60 

+72 53 43.2 

8.01 

7 

1 

G110.80+14.16 

400 

14+1 

21 57 35.84 

+72 46 49.1 

6.88 

2 

1 

G111.41-2.95 

3000 

13+2 

23 22 28.78 

+57 36 44.5 

7.77 

1 

1 

G131.65+9.75 N 

1070 

11+1 

02 40 03.30 

+70 43 20.1 

6.79 

4 

1 

G131.65+9.75 S 

1070 

13+1 

02 40 11.20 

+70 36 09.5 

8.89 

7 

1 

G132.12+8.95 

850 

10.6+0.9 

02 41 49.64 

+69 50 12.9 

5.92 

2 

1 

G149.67+3.56 

170 

13+1 

04 17 08.20 

+55 13 40.6 

4.89 

1 

1 

G154.08+5.23 

170 

11.2+0.9 

04 48 08.73 

+53 07 23.1 

20.2 

3 

1 

G157.92-2.28 

2500 

10.9+0.8 

04 28 48.30 

+45 24 22.5 

8.33 

2 

1 

G159.34+11.21 

1590 

14+1 

05 42 17.30 

+52 08 16.4 

1 point 

1 

1 

G161.55-9.30 

250 

13+1 

04 16 15.30 

+37 45 35.0 

2.22 

5 

1 

G164.71-5.64 

330 

13+1 

04 40 44.43 

+37 45 14.3 

5.39 

1 

1 

G167.20-8.69 

160 

12+1 

04 35 26.88 

+34 18 25.6 

6.88 

2 

1 

G168.85-10.19 

2610 

13+1 

04 37 07.59 

+31 43 01.2 

1 point 

1 

0 

G173.43-5.44 

150 

13+1 

05 07 49.46 

+31 20 46.3 

1 point 

1 

0 


2.3. Dust continuum observations with Herschel 


Table 2. Center coordinates of the 1.5' radius reference areas. 


We observed the dust continuum data with the Herschel SPIRE 
instrument ( [Griffin et al.||2010} > between November 2009 and 
May 2011. These observations consist of surface brightness 
maps of 250, 350, and 500 (rm. The raw and pipeline-reduced 
data are available via the Herschel Science Archive, the user- 
reduced maps are available via ESA sit^j 

We reduced these data similar to ( [Juvela et al.||2012b| and 
only a summary of these reduction steps is given here. We com¬ 
pleted the reduction using the Herschel Interactive Processing 
Environment HIPE v.10.0 and the official pipeline with the it¬ 
erative destriper and extended emission calibration options. The 
resulting maps are the product of direct projection onto the sky 
and averaging of the time-ordered data. The resolutions of the 
maps are 18", 25", and 37" at 250, 350, and 500 inn, respec¬ 
tively. 

The accuracy of the gain calibration of Herschel data is bet¬ 
ter than 7 % in absolute terms and probably better than 2 % band- 
to-band Q The surface brightness data do not have an absolute 
zero point and therefore, before temperatures or column densi¬ 
ties can be estimated, we must set a consistent zero point across 
all Herschel bands. One alternative is to compare Herschel data 
with Planck and IRAS measurements (using IRAS data tied to 
DIRBE scale) and to make use of the zero points of those surveys 
(see, e.g., Juvela et al .|2012b) . We opt for the more straightfor¬ 
ward procedure of selecting a reference area with a radius of 1.5' 
within the Herschel map and measuring surface brightness val¬ 
ues relative to the average value found in the reference area. The 
central coordinates of the reference areas can be seen in Table 
[2] The derived dust temperature and column density estimates 
ignore the very diffuse medium to the extent that is visible in 
the reference region. The reference regions have low extinction 
(Ay ~ 2 or less) and therefore probably a lower fractional abun- 


2 http://herschel. esac. esa. int/UserReducedData. shtml 

3 SPIRE Observer’s manual 
http://herschel. esac. esa. int/Documentation. shtml 


Field 

a (J2000) 

8 (J2000) 

G86.97-4.06 

21 16 40.8 

+43 34 45 

G92.04+3.93 

21 01 17.9 

+52 34 45 

G93.21+9.55 

20 37 44.3 

+56 44 49 

G94.15+6.50 

20 59 48.7 

+55 56 19 

G98.00+8.75 

21 05 59.4 

+59 51 19 

G105.57+10.39 

21 43 14.5 

+66 26 12 

G108.28+16.68 

21 08 16.0 

+72 51 20 

G110.80+14.16 

21 59 05.7 

+72 37 40 

G111.41-2.95 

23 20 11.9 

+57 35 20 

G131.65+9.75 

02 42 10.5 

+70 45 20 

G132.12+8.95 

02 37 04.7 

+69 53 43 

G149.67+3.56 

04 18 42.2 

+55 29 55 

G154.08+5.23 

04 47 58.8 

+53 14 42 

G157.92-2.28 

04 30 02.5 

+45 25 50 

G159.34+11.21 

05 41 07.7 

+51 54 33 

G161.55-9.30 

04 15 46.2 

+37 40 56 

G164.71-5.64 

04 42 47.6 

+38 22 45 

G167.20-8.69 

04 35 10.1 

+34 11 09 

G168.85-10.19 

04 36 06.8 

+31 50 54 

G173.43-5.44 

05 07 55.2 

+31 04 55 


dance of CO. Thus, even after background subtraction, the con¬ 
tinuum data may probe a volume larger than the actual molecular 
cloud. Therefore, the extrapolation of the relations like in Fig. |T] 
does not necessarily go via the origin of that plot. This should 
not be a significant problem, however, because we are interested 
in the densest clumps that are far above the column density of 
the reference regions. 

The resolution of Herschel column density maps is 40". The 
beam size is 35" in CO observations and 40" in N 2 IF obser¬ 
vations. Because of the relatively small differences and the ex¬ 
tended nature of our sources, the data were compared directly 
without further convolution. 
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3. Methods 

For the calculation of the CO column densities, we used the 
method described by Myers et al. ( | 1983) >. Based on the data anal¬ 
ysis, we calculated the optical depth at the peak of C 18 0 line, ns, 
the excitation temperature of the C 18 0 line, 7j x , and the column 
density of the C 18 0 in the cloud, /Vis. To derive the total column 
density of Fh, the column density of C 18 0, Ni%, is divided by the 
commonly assumed abundance ratio [Cl 1 s O]/| H 2 ] = 10 1 . 

The optical depths of 13 CO and C 18 0 are related by (Myers 
etaLl l fl9gi3) , Eg. 2) 

n 13 (J = 1) L 13 av 18 /(r 18 ) 

Tn = Tio-. (1) 

nit(J = 1) L W AV 13 J(T u ) 

In the equation J(T) = 7o[exp(7o/7 ) - 1] _1 , where 7 0 = 5.27 K 
for C 18 0 and 5.29 K for 13 CO and T is the excitation temper¬ 
ature. The 7 ?i 3 (J=l) and 77 is(J=l) are the number of molecules 
at J=1 levels of the 13 CO J=l-0 and C 18 0, respectively. The Li 3 
and Lis are the line-of-sight extent of the emitting gas and the 
AV 13 and A Fix are the line widths of the 13 CO and C ls O line, 
respectively. The equation assumes the same excitation temper¬ 
ature for both lines, the same beam filling (both lines originate 
in the same region), and the same velocity gradient. We also as¬ 
sume the terrestrial abundance ratio of 5.5 between l3 CO and 
C 18 0 and we use the line width of C 18 0 for the calculations. 


For a more comprehensive explanation of the method, see Myers 
|etaL[ ( fT983| . 

We also derived the column density with the Herschel dust 
continuum observations. The observed intensity can be stated as 


I v = BvCT’du.tXl -e T ) « B(7’ dust )xr, 


( 2 ) 


T = KyX ATH 2 ) X fl. 


(3) 


where 11 is the average particle mass per IT molecule, 2.8 u. 
We assume a dust opacity of k = 0.1(//1000.0 x 10 9 ) 20 cm 2 /g, 
where / is the frequency ( [Beckwith et al.||T990 >. The constant 
2.0 is the dust opacity spectral index, /7. While it may vary 
from source to source, the value 2.0 we used is consistent with 
many observations of dense clumps although the average value 
in molecular clouds is likely to be closer to ~1.8 ([Beckwith 
et al.||1990[ |Boulanger et al ||1996[ |Hildebrand|fl983[ |Planck 
Collaboration et al.|201 la|c|e) . Again, the column density of hy- 
drogen is A(H 2 ). From Eqs.|2|and[3]we can derive the molecular 
column density 


ATH 2 ) = 


(4) 


77,(7 t ii ls t )K..p 

The derived column densities from molecular lines and dust 
can be used to calculate the mass and density of a cloud or a 
clump as 


M c = A(Hi) x nR 2 x —, 
M g 


(5) 


where 7? is the radius of the clump. The average density, n, is 
calculated using the hydrogen column density 


MH 2 ) 
27? ’ 


where the clump diameter 27? is the full-width at half maximum 
(FWHM) from the Herschel column density map. 

We derived the estimated virial mass from the formula 
( jMacLaren et al.|1988) > 


kef R 


(7) 


where k depends on the density distribution. We choose k = 
1.333, which corresponds to density distribution p(r) oc 5 . 
The velocity dispersion is given by the equation 


\kTkin | 

I AV 2 kT kin \ 

I m 

(8/77 ( 2 ) 777 ) 


( 8 ) 


where m is the mean molecular mass (2.33 u assuming 10 % 
He), and m is the mass of the molecule used for observations. 
The temperature, T kln , is the kinetic temperature, which is as¬ 
sumed to be 10 K. The velocity dispersion provides an estimate 
of the nonthermal motions that provide further support against 
gravity (Planck Collaboration et al. [2011b , Bertoldi & McKee 
|1992| ). If the cloud’s mass is less than the virial mass, it is not 
gravitationally bound, and without external pressure for support, 
it will disperse. 

We can also examine the stability using the model of BE 
spheres. The BE mass assumes a static isothermal cloud with no 
magnetic field ( McKee & Qstriker||2007 1 . We include the non¬ 
thermal component and use the effective sound speed instead of 
the isothermal sound speed. As a difference to the virial mass, the 
BE model includes the external pressure. The critical BE mass 
dBonnor|1956| ) is 


where I v is the observed intensity at a frequency v, 7? v (7’ dust ) is 
the blackbody brightness of the object as a function of color tem¬ 
perature Tdust , and r is the source optical depth. The equation 
assumes a homogeneous source. In the far-infrared and at longer 
wavelengths the optical depths of the clouds are clearly below 
one. This justifies the approximation made in Eq.[2] The optical 
depth can be written as 


M, 


C7~ 


BE 


2.47? BE —, 


(9) 


(6) 


where 7 ?be is the BE radius, and G is the gravitational constant. 

4. Results 

4.1. Comparison of dust continuum and moiecuiar tine data 

First we compare the column densities calculated from the 
molecular line observations to the column densities calculated 
from the dust continuum observations. The sizes of the clumps 
were estimated as a FWHM from the dust column density maps 
and the physical sizes in pc are in Table [3] Twelve of the clumps 
are larger than the size criteria (0.03 - 0.2 pc) given by Bergin 
[& Tafalla] |2007| for cores, and, thus, are larger structures that 
could contain further substructures. 

We did not get a result for the excitation temperature in all of 
the cases. This is most likely because of the greater ratio between 
I3 CO and C 18 0 abundances than the 5.5 that was assumed in the 
calculations. Thus, we initially used an estimate of T ex =1 0 K 
(see, e.g., |Dobashi et al7||1994 |Alves et al.||l999| . The results 
from the calculations showed, however, that when T ex could be 
estimated, it was typically 5 K (see Table 3]). Therefore, we show 
in Table[3]column densities for a fixed value of T ex =5 K. We plot 
the column densities calculated from molecular observations in 
comparison with the column density calculated from dust con¬ 
tinuum observations in Fig. |T| 

In Fig. |T] the cases where 7' cx could be derived are plot¬ 
ted with square symbols, the solid line showing the corre¬ 
sponding linear least-squares fit. The column densities calcu¬ 
lated with T ex =5 K are plotted with circles and fitted with 
a dashed line. In the cases where the excitation tempera¬ 
ture could be calculated from the observations, the T ex -5 K 
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Table 3. Clump sizes (FWHM from dust column density), column densities derived from dust and line data, and estimated C ls O 
excitation temperatures. 


Field 

FWHM [’] 

FWHM [pc] 

m h 2 ) 

dust [10 21 cm 2 ] 

Tex [K] 

N( H 2 ) 

lines [10 21 cm 2 ] 

N( H 2 ), r cx =5 K 
lines [10 21 cm 2 ] 

G86.97-4.06 

0.4 

0.17 

12±4 

- 

- 

4+1 

G92.04+3.93 

0.8 

0.38 

50+20 

4+2 

90+40 

37+8 

G93.21+9.55 

0.6 

0.15 

19+6 

4.9+0.4 

9+2 

8+1 

G94.15+6.50 

0.9 

0.41 

7+2 

4+1 

19+8 

13+2 

G98.00+8.75 

0.5 

0.35 

10+4 

6+1 

6+1 

6.5+0.9 

G105.57+10.39 

0.5 

0.28 

15+5 

- 

- 

8+1 

G108.28+16.68 

0.8 

0.15 

1.5+0.4 

- 

- 

2.7+0.8 

G110.80+14.16 

1.0 

0.24 

2.2+0.7 

4.7+0.3 

3.3+0.8 

3.0+0.3 

Gil 1.41-2.95 

0.6 

1.09 

6+2 

- 

- 

4.2+0.1 

G131.65+9.75 N 

0.6 

0.38 

9+3 

- 

- 

5+1 

G131.65+9.75 S 

0.3 

0.17 

4+1 

- 

- 

3.6+0.5 

G132.12+8.95 

1.3 

0.64 

15+5 

- 

- 

5+1 

G149.67+3.56 

0.4 

0.04 

5+2 

- 

- 

1.6+0.4 

G154.08+5.23 

0.9 

0.09 

17+6 

5.5+0.5 

9+2 

11+1 

G157.92-2.28 

0.4 

0.55 

5+2 

- 

- 

1.8+0.6 

G159.34+11.21 

1.4 

1.33 

4+1 

- 

- 

6+2 

G161.55-9.30 

0.5 

0.07 

6+2 

- 

- 

1.6+0.4 

G164.71-5.64 

0.8 

0.15 

3+1 

- 

- 

1.7+0.3 

G167.20-8.69 

1.2 

0.12 

5+2 

5.1+0.4 

4+1 

4.7+0.5 

G168.85-10.19 

1.6 

2.43 

1.3+0.4 

- 

- 

0.2+0.1 

G173.43-5.44 

2.6 

0.23 

3+1 

- 

- 

1 . 0 + 0.2 



Fig. 1 . Comparison of the column densities derived from molec¬ 
ular line and dust continuum observations. The squares are the 
column densities, where the excitation temperature could be 
calculated, and the circles are the column densities calculated 
with an assumed excitation temperature of 5 K. The colors in¬ 
dicate the dust color temperature. Blue is below 11 K, yellow 
is 11 - 13 K, and red is warmer than 13 K. The solid line 
is the linear least-squares fit to the column densities derived 
from the calculated excitation temperature, and the dashed line 
fits the T ex -5 K values. The green line indicates a line where 
JV(H 2 ) dust =ATH 2 ) lme . 


values are usually within the margin of errors. The col¬ 
umn densities calculated from dust are typically higher but 
also often within the margin of error. For example, for 
G159.34+11.21, for whiclT A(H 2 ) dust =4±lxl0 21 cnr 2 and 
^V(H 2 )5 k=6±2x 10 21 cm -2 . In some cases the column densities 
are quite different, however. This is the case for G86.97-4.06, 


where ATH 2 ) duM = 1,2±0.4x 10 22 cm -2 , while the value derived 
from the line data is much lower, ATH 2 )5 k= 4±1 xlO 21 cm~ 2 . 
Uncertainty of the background subtraction could contribute to 
some of these differences. For example, in G86.97-4.06 the 
reference region is estimated to have a column density of 
~10 21 cm 2 , which corresponds to a diffuse part of the cloud with 
little CO. However, the presence of a small column density of 
gas without CO molecules cannot explain the difference in the 
column density estimates of several times 10 21 cm 2 . 

The dust color temperatures of the clumps are shown in Fig. 
[T]in different colors: blue for clumps below 11 K, yellow for the 
range of 11-13 K, and red above 13 K. The column densities 
derived from dust show an anticorrelation with the dust color 
temperature as denser regions are cooler. This anticorrelation on 
the dust color temperature was not seen as clearly in the column 
densities derived from molecular line observations, as there are 
several column densities ~ 0.5 x 10 22 cm -2 independent of the 
dust color temperature. 

The maps of 7 m b( l , CO) are shown in Figs. [2] and [3] The 
contours correspond to the column densities derived from dust 
continuum data (subtracting the local diffuse background). The 
statistical error in the column density maps of the dust is at level 
2xl0 20 cm~ 2 . All the contours are at more than 10cr level and, 
thus, statistically significant. The S/N of l3 CO maps is lower. In 
a few cases, we get over 5<x detection for the whole map, namely 
for fields G105.57+10.39 and G149.67+3.56. In five fields, 
G86.97-4.06, G108.28+16.68, G131.65+9.75 S, G157.92.28, 
and G164.71-5.64, the signal drops below 2cr in part of the field. 

In all the maps, we see a peak in the main beam temperature 
of 13 CO close to the dust column density peak. For few fields, 
e.g., G98.00+8.75 (Fig.[2]», the peaks coincide, but for most cases 
there is a small shift, which can be up to 0.5'. In most cases 
the shift is smaller than the resolution of the observations. In 
some cases, as in G149.67+3.56, the peaks coincide partly, but 
are still shifted. The shifts are small in physical sizes as well, 
~0.01 -0.22 pc, and less than the radius of the clumps. 
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10.00s lOmOO.OOs 50.00s 21h09m40.00s 

RA (J2000) 




30.00s 27.00s 24.00s 21.00s 21h02ml8.00s 

RA (J2000) 









Fig. 2. Maps of l3 CO main beam temperature (r m b [K], resolution 35")- The contours show the column density derived from dust 
continuum observations (resolution 40"). with contour steps 10 % of the peak value. The name of the field and the maximum value 
of N(H 2 ) derived from dust emission are marked in each frame. 


Significant differences in the peak positions of dust and line 
emission are often observed because of molecular depletion, 
other chemical gradients, or because of temperature differences. 
[Marsh et al.[ ( |2014[ > found shifts with a median distance of 59" 
for H 3 CO 1 and dust (250 (im). The dust and gas have been 
found to be coupled in the interior of cores with densities above 
3xl0 4 cm -3 . At lower densities the temperature difference be¬ 


tween the two components depends on the strength of the ex¬ 
ternal UV field and not on density ( jGalli et al.[2002) l. Also, the 
depletion can increase the gas temperature in low and moderate 
densities ( Goldsmith|2001| >. If we compare the densities calcu¬ 
lated from the molecular line observations in Table[4j the clumps 
are less dense and, thus, dust and gas are not likely coupled 
in these clumps. However, the asymmetric displacement of the 
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Fig. 3. Continued from Fig. [ 2 ] 


l3 CO emission in the densest clumps is more likely to be a sign 
of variations of 13 CO abundance. 


4.2. Clump masses 


The masses and densities along with the virial and BE masses 
are listed in Table [4] The uncertainties for densities and masses 
were determined with the Monte Carlo method. The uncertain¬ 
ties are big in many cases, up to 53 %, and we must be care¬ 
ful when drawing conclusions from the density and masses. The 
mass errors do not include the error for distance (~30 %, in the 
case of G168.85-10.19 the error is closer to ~50 %), which fur¬ 
ther increases the uncertainty. In addition, in most cases the ex¬ 
citation temperature could not be calculated and the excitation 
temperature of 5 K was assumed. In the cases where the excita¬ 
tion temperatures could be calculated, the derived densities and 
masses are within the margin of error of those derived from the 
dust emission data. The field G92.04+3.93 is the only exception. 

The calculated masses vary greatly between the clumps, 
from one tenth of a solar mass to a couple of hundred. Thus, 
they clearly cover a mix of different structural scales. Following 
the mass criteria of Bergin & Tafalla| ( f2007] l, 11 of them are more 
massive than typical cores (> ocal5M Q ). However, most of them 
show a single kinematic component and are therefore consistent 
with the idea of clumps as velocity coherent objects. If both Mi; ne 
and M dust are larger than the virial and BE mass, we consider 
them to be gravitationally bound, and thus (potentially) on the 
verge of collapse. If M| mc or M dus , exceeds BE mass, the object 
is potentially bound, but if the estimated mass is smaller than 
both the virial and BE mass, we consider the object not to be 
prestellar. 

Two clumps, G92.04+3.93 and G159.34+11.21, appear to 
be on the verge of collapse. The mass of G92.04+3.93 calcu¬ 
lated with T ex =5 K is lower than the virial mass, but when exci¬ 
tation temperature could be calculated it exceeds the virial mass 


clearly. For the held G159.34+11.21, the excitation temperature 
could not be calculated, thus, the calculated mass is based on 
the assumption of T ex =5 K. Three other clumps, G94.15+6.50, 
G98.00+8.75, and G105.57+10.39, are supercritical, albeit not 
significantly, if we consider the BE mass. The remaining 16 
clumps seem to be near equilibrium. 


4.3. Observations ofN 2 H + 


To estimate the NTH 1 column densities, we performed fits of the 
hyperhne structure of the J =1-0 lines. The fits were carried out 
assuming a single excitation temperature. We treated the excita¬ 
tion temperature, optical depth of the 23-12 component, and the 
radial velocity and width of the spectral lines as free parame¬ 
ters (see Pirogov et al. 2003[ >. Because of the low signal-to-noise 
ratio, the optical depths and excitation temperatures were un¬ 
certain. Therefore we calculated the column density assuming a 
fixed value of the excitation temperature, T ex =5 K, and optically 
thin emission (Caselli et al. 2002) . The line parameters and ob¬ 
tained column densities are listed in Table 0 We estimated the 
error estimates with the Monte Carlo method, however, exclud¬ 
ing the effect of the uncertainty of T ex . Apart from G86.97-4.06, 
shown in Fig. [4] the N 2 H 1 spectra are shown in Appendix [b] 

For most of the studied clumps, N 2 H + line was not de¬ 
tected, suggesting that these clumps are not dense enough for 
significant CO depletion and the associated rise of N 2 H 1 abun¬ 
dance (|Bergin & Langer|l997||Charnley|l997|). The line was de¬ 
tected in six clumps: G86.97-4.06, G92.04+3.93, G93.21+9.55, 
G98.00+8.75, G105.57+10.39, and G132.12+8.95 (see Fig. |T] 
for an example of the spectra and Table [5j. These clumps were 
among those with coldest dust temperatures. The main line in 
the fields was T m b ~0.2—0.8 K and the satellites are barely visi¬ 
ble above the noise of -0.05 K. The low excitation temperature 
of the CO lines can also be an indication of a low density, as in 
lower densities the molecules are less thermalized and the exci- 
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Table 4. The used line width of C 18 0, densities, masses, virial masses, and BE masses for all the targets using molecular line 
observations and mass value calculated from the dust observations. The values from molecular line observations are derived with 
T ex calculated from C 18 0 observations and with a fixed value of T ex = 5 K. 


Field 

AVpiso 

[km s -1 ] 

71 [10 3 cm- 3 ] 

n [10 s cm 3 ] 
(Tex-5 K) 

Mune [M 0 ] 

Mu™, [M 0 ] 
(Tex=5 K) 

M dust [M 0 ] 

M vir [M 0 ] 

M be [M 0 ] 

G86.97-4.06 

0.74 

- 

7±2 

- 

1.7+0.5 

6+2 

10+2 

6+2 

G92.04+3.93 

1.79 

70±20 

32±8 

230+70 

100+20 

130+40 

110+30 

70+10 

G93.21+9.55 

0.65 

19±4 

18±2 

3.2+0.7 

3.1+0.4 

7+2 

7.5+0.8 

4.5+0.5 

G94.15+6.50 

0.85 

15±8 

11±2 

60+30 

39+7 

20+7 

31+6 

18+3 

G98.00+8.75 

0.76 

5+1 

6.0+0.8 

12+2 

13+2 

22+9 

22+3 

13+2 

G105.57+10.39 

0.76 

- 

9±1 

- 

11+2 

21+7 

18+3 

11+1 

G108.28+16.68 

1.07 

- 

6±2 

- 

1.0+0.3 

0.6+0.2 

16+4 

10+2 

G110.80+14.16 

0.31 

4.4+0.9 

4.0+0.5 

3.3+0.7 

3.0+0.4 

2.2+0.7 

5.6+0.6 

3.4+0.4 

Gil 1.41-2.95 

2.27 

- 

1.2+0.3 

- 

90+20 

120+50 

500+100 

290+80 

G131.65+9.75N 

1.22 

- 

4±1 

- 

13+3 

23+8 

50+10 

30+8 

G131.65+9.75 S 

0.53 

- 

7±1 

- 

1.7+0.2 

1.8+0.6 

6.5+0.9 

3.9+0.6 

G132.12+8.95 

1.39 

- 

2.3+0.5 

- 

32+7 

110+30 

110+30 

70+20 

G149.67+3.56 

0.64 

- 

12±3 

- 

0.05+0.01 

0.15+0.07 

2.1+0.5 

1.3+0.3 

G154.08+5.23 

0.56 

33±7 

40±4 

1.3+0.3 

1.5+0.2 

2.3+0.8 

3.7+0.5 

2.2+0.3 

G157.92-2.28 

0.87 

- 

1.1 ±0.3 

- 

9+3 

30+10 

40+10 

26+8 

G159.34+11.21 

0.84 

- 

1.4+0.4 

- 

180+60 

130+30 

100+30 

60+20 

G161.55-9.30 

0.48 

- 

7±2 

- 

0.13+0.04 

0.4+0.1 

2.3+0.7 

1.4+0.4 

G164.71-5.64 

0.50 

- 

3.6+0.7 

- 

0.7+0.1 

1.4+0.4 

6+1 

3.4+0.7 

G167.20-8.69 

0.32 

12±2 

12±1 

1.0+0.2 

1.1+0.1 

0.47+0.05 

2.8+0.3 

1.7+0.2 

G168.85-10.19 

0.21 

- 

0.03+0.01 

- 

20+10 

130+40 

50+30 

30+10 

G173.43-5.44 

0.24 

- 

1.4+0.2 

- 

0.9+0.1 

2+1 

4.5+0.7 

2.7+0.4 


Table 5. The main beam temperature, velocity, noise and column density of N 2 H + in the fields it was detected. 


Field 

T mb [K] 

V [km s 1 ] 

AV [kms- 1 ] 

rms [K] 

N [10 12 cm--] 

G86.97-4.06 

0.34+0.03 

5.50+0.02 

0.43+0.04 

0.06 

1.1+0.3 

G92.04+3.93 

0.17+0.01 

-1.80+0.02 

0.51+0.04 

0.04 

3.0+0.2 

G93.21+9.55 

0.23+0.04 

-1.73+0.03 

0.41+0.08 

0.07 

2.1+0.5 

G98.00+8.75 

0.23+0.02 

4.98+0.03 

0.54+0.04 

0.06 

2.3+0.6 

G105.57+10.39 

0.31+0.02 

-10.10+0.02 

0.46+0.04 

0.05 

2.0+0.5 

G132.12+8.95 

0.32+0.02 

-12.41+0.02 

0.51+0.04 

0.04 

1.1+0.3 


tation temperature is much lower than the kinetic temperature. 
In the studied clumps, the excitation temperature of the C ls O in 
most cases was ~5 K indicating a clear difference between the 
excitation and kinetic temperatures. 


The calculated N 2 I-E column densities are only a few times 
10 12 cm -2 . These can be compared with the N 2 H + survey of 
nearby, low-mass cores conducted by Caselli et al.]( 2002 1 . In that 
study, for the low-mass cores, M vn - ~ 10 Mq, the average N 2 H + 
column density was ~7x 10 12 cm -2 , and they found starless cores 
had a marginally lower column density than protostellar cores. 
Our column densities correspond to a lower end of those val¬ 
ues. Conversely, Pirogov et al. (2003) observed a sample of 35 
massive molecular cloud cores (mean virial mass over 600 solar 
masses) and found an average column density of 29xl0 12 cm -2 . 
Our low column density values are consistent with the low mass 
of the clumps, but also suggest that the densities (and CO deple¬ 
tion) are not significant enough for a large N 2 PC abundance. 


4.4. Modeling 

The observations were compared to radiative transfer models of 
CO lines and of dust continuum emission. The modeling is ex¬ 
plained in more detail in Appendix [C] and we present only the 
main results. The modeling was limited to the two positions in 
the field G131.65+9.75. These positions were chosen as they 


were in the same field and, thus, enable comparison of clumps at 
the cloud edge and in its center. 

The calculations that are based on the dust model of 
Ossenkopf & Henning (1994) find 13 CO abundances that are 


close to the value of 10 - °. Because the abundances depend 
greatly on the column density estimated from continuum data, 
the values are uncertain. The abundance ratio for 13 CO and C 18 0 
was found to be ~10, instead of the 5.5 assumed in Sect. [3] 
independent of the different assumptions used in the modeling 
(see Table |C.1~| . Harjunpaa & Mattila (1996) and Anderson et al.| 


have also found similar abundance ratios in dark clouds 
and cores. 


5. Discussion 

In about one third of the clumps, where excitation temperature 
could be estimated, the column densities are similar to those de¬ 
rived from dust (see Fig. [Tjand Table [3]). The column densities 
calculated with r ex =10 K, were lower by a factor of three. Using 
T ex =5 K, which was more in line with those that could be cal¬ 
culated, we get very similar results, especially with the cases 
where column densities are lower (~10 21 cm -2 ). Higher column 
densities tend to have larger uncertainties. For example, the col¬ 
umn density for clump G92.04+3.93, 5±2xl0 22 cm -2 from dust 
and 9+4xl0 22 cm 2 from molecular lines is still inside the mar¬ 
gin of errors. The cases where T ex =5 K were assumed to have 
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Fig. 4. The N 2 H + spectrum observed toward G86.97-4.06. The 
red dashed line is the fit to the hyperfine spectra and the lower 
subframe shows the fit residuals. The main beam temperature 
of the 23-12 component, the fitted radial velocity, the FWHM 
line width (km s 1 ), and the residual rms noise (as main beam 
temperature) are given in the frames. 



Fig. 5. FWHM (pc) vs. distance with the best fit. 


mostly lower column densities than those where T ex was solved. 
One explanation for this could be that CO has been depleted in 
these areas and, thus, the column density and the volume density 
calculated from CO observations are too low. The most likely 
explanation for the difference to the dust estimates is, however, 
the excitation temperature that is likely to vary from clump to 
clump. 

The distances vary from 150 pc to 3 kpc across the sample, 
but the physical sizes (see Table [3] and Fig. [5| of the objects 
are, with one exception, ~ 1 pc or below and thus comparable to 
typical clump scales as defined by jBergin & Tafalla| ( |2007| ). The 
largest sizes include those of the most distant clouds, namely. 
Gill.80+14.16, G157.92-2.28, G159.34+11.21, and G168.85- 
10.19. On average, the physical size is proportional to 0.77 times 


the distance, a natural result of our resolution at larger distances 
corresponding to larger structures. 

Although there are large uncertainties in the column densi¬ 
ties derived from the lines, even with the calculated T ex values, 
the mass and column density estimates indicate that the selected 
clumps might not be very dense, and subsequently may not be 
actual prestellar cores. The lack of N 2 H + emission also supports 
this conclusion, with upper limits 7 mh (N 2 H 1 )=0.05—0.2 K for 
the 16 of 21 clumps. 

The fact that the excitation temperature could not be calcu¬ 
lated in all of the selected positions could be caused by the abun¬ 
dance ratio of l3 CO and C 18 0 , which was actually higher than 
the assumed upper limit of 5.5. When extinction Ay is small, 
the abundance ratio can range much above the terrestrial value 
dMinchin et al.||1995) >. The abundance ratio we found in our 
models was ~10, compared to the used terrestrial value of 5.5. 
We found the abundance for l3 CO, however, to be close to the 
canonical 10 6 , provided that we trust the dust-derived column 
densities. The difference in the estimated abundance ratio can 
be caused by a real differences in the abundances or excitation 
temperatures or by optical depth effects. The modeling, how¬ 
ever, roughly takes the effect of optical depth and the expected 
difference in the excitation temperature of the isotopomers into 
account. This difference is caused by the difference in the pho¬ 
ton trapping resulting from the difference in the optical depth of 
the lines. The assumption of the homogeneous source in the LTE 
analysis is inaccurate and, in reality, the excitation temperature 
of C ls O is probably below the excitation temperature of l3 CO 
(and not equal as assumed) and drops precipitously toward the 
cloud surface. If the 13 CO is optically thick, the observed in¬ 
tensity originates in a different part of the cloud than the C 18 0. 
Our calculations did not show 13 CO to be optically thick, but 
we calculated the optical thickness using the assumption that the 
excitation temperatures are equal. Because of the smaller optical 
depth, more of the C 18 0 line photons escape, leading to a smaller 
excitation temperature. This is also clearly significant when the 
kinetic temperature varies within the cloud (Juvela et al .|2012aj >. 

Based on the line data, 16 of the studied clumps are sub- 
critical and, even considering the uncertainties of the mass esti¬ 
mates, only five sources could be gravitationally bound. In three 
cases, G92.04+3.93, G94.15+6.50, and G159.34+11.21, the cal¬ 
culated mass is larger than the virial and BE mass even when 
considering the uncertainty. Two other clumps, G98.00+8.75 
and G105.57+10.39, are within error margins when we con¬ 
sider the BE mass. The five clumps are from the entire temper¬ 
ature range. However, only one clump, namely G92.04+3.93, 
had a mass larger than the virial and BE masses on a level 
> lcr, although G92.04+3.93 does not exceed this level if 
we assume 7' ex =5 K. If we consider the error in the dis¬ 
tance of the object, none of the sources exceed the lcr-level. 
However, only four clumps have masses below the virial and 
BE masses on a level > lcr. In the coldest supercritical clumps, 
G92.04+3.93 (T dust =10.4 K), G98.00+8.75 (r dust =11.6 K), 
and G105.57+10.39 (T dust =10.9 K), we detect the N 2 H + line, 
which would support the idea that these are being denser and 
possibly prestellar. The N 2 H + emission seemed in general to 
be connected to the coldest clumps of —10— 11 K, although it 
was not observable in all the clumps of these temperatures. For 
the clumps where N d H + was detected, we get a FWHM line 
width of <0.5 km s 1 , values typically found in prestellar cores 
(Johnstone et al.|2010| . 

If we rely on the mass estimates derived from dust ob¬ 
servations instead of CO lines, the results are similar, only 
three clumps are below the virial mass and BE mass limits. 
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Four of the clumps were larger by more than the error mar¬ 
gins and could be classified as prestellar. The four clumps were 
G92.04+3.93, G105.57+10.39, G159.34+11.21, and G168.85- 
10.19. Thus, most of the clumps that were detected for low dust 
color temperature appear to disperse with time rather than form 
. Meng et al. (2013| and Wu et al. (|2012| had similar results 


stars. 


in their cold cores surveys: only a small fraction of the clumps 
were found to be prestellar. Of the 673 sources, 10 clumps were 
mapped and 22 potential cold cores were identified by Wu et al. 


( |2012| >. Of these 22 cores, seven were found to be gravitationally 
bound. Our results are in general agreement with these results 
on the fractions of starless and prestellar cores (5 out of 21 were 
gravitationally bound). Meng et al. ( 201 3| > studied a subset of 
|Wu et al.] ( |2012| l in more detail, and 90 % of them were found to 
be starless (4 cores out of a sample of 38 were associated with 
sources). Studies of protostellar sources should be extended to 
a larger set of Planck clumps, including the 21 clumps that we 
discuss. 


Appendix B: N 2 H + 

The spectra with N 2 H + detections, besides the field G86.97-4.06 
in Fig. [4] are shown in Fig.|B.l| 


Appendix C: Radiative transfer models 

In Sect. [4] the dust emission was analyzed assuming a constant 
temperature along the line-of-sight (LOS), but this may under¬ 
estimate the true column density of externally heated, optically 
thick clumps ( Malinen et al.|20lT| Juvela et al.|2013| l. Similarly, 
the 13 CO(1-0) and C ls O(l-0) lines were analyzed using the LTE 
assumption, without the possibility of independent estimates of 
their relative abundances. To examine these questions, we car¬ 
ried out radiative transfer modeling where the density distribu¬ 
tion was first derived from dust continuum observations and was 
then used as a basis for modeling the lines. The modeling was 
limited to the two positions in the field G131.65+9.75. 


6. Conclusions and summary 

We investigated 21 clumps, that show low dust color tempera¬ 
tures T~ 10-15 K and are therefore potential places of star for¬ 
mation. To study the conditions in the clouds, we used molecu¬ 
lar line and dust continuum observations. Our comparison of the 
two tracers shows that, even though the gas and dust emission are 
mainly morphologically compatible, the dust peak is sometimes 
shifted up to 30" relative to the 13 CO maximum. The data were 
examined with standard LTE analysis and with radiative transfer 
modeling to understand better the physical and chemical state of 
the clumps. 

The column density calculations from molecular lines were 
mostly within the uncertainties to those calculated from dust, 
when the excitation temperature could be calculated. With an as¬ 
sumed fixed excitation temperature value of T ex =5 K, the derived 
column densities were only slightly below the dust estimates. 

Using modeling to compare dust continuum and line data, 
the abundance of 13 CO was found to be close to the typically 
assumed value of 10 6 . However, the abundance ratio of l3 CO 
and C 18 0 was ~10, higher than the terrestrial value 5.5. 

The calculated masses had big uncertainties and one must 
be careful when drawing conclusions. When the masses were 
compared to virial and BE masses, only five clumps had a 
high enough mass to be gravitationally bound. In three of these 
clumps, we also found N 2 H + , so the clumps found in the 
fields G92.04+3.93, G98.00+8.75, and G105.57+10.39 could be 
prestellar. G92.04+3.93 and G105.57+10.3 were also among the 
coldest T dust ~ 10-11 K, while G98.00+8.75 was slightly warmer 
Vdus, =12 K. The highest temperatures in the examined clumps 
were 7d ust >13 K. 

The stability of the clumps requires further study. 
Observations of higher CO isotopomer transitions and further 
density and temperature tracers are needed to better quantify the 
gas component in these cold clumps. 


Appendix A: 13 CO and C ls O spectra 


isr 


A sa mpl e spe ctrum of l3 CO of all the fields are shown in Fig. 


A.l 


and 


A.2 


The spectra are from the l3 CO peak position of 
each field. The C 18 0 spectrum is shown from the same position; 
the baseline has been moved to -2.0 K. 


C. 1. Modeling of dust surface brightness 


We extracted 5'x5' continuum surface brightness maps centered 
on the selected positions and resampled the data on 3" pixels. 
The maps are 100x100 pixels in size and corresponding model 
clouds were constructed using a cartesian grid of 100 x 100 cells. 
In the line-of-sight direction the density distribution was as¬ 


sumed to follow a Plummer-like (Whitworth & Ward-Thompson 


2001 


Plummer 


1911 


fiat part with K 


function p(r) 


1 r 

[|l+(r/*na.) 2 |] 


with a central 


(i at equal to 0.03 pc (not well resolved with the 
employed discretization). Corresponding to the apparent clump 
sizes in the plane of the sky, the FWHM of the density distri¬ 
bution was set to ~17 pixels, which corresponds to ~ 50" or a 
linear scale of 0.25 pc at the distance of 1070 pc. This may over¬ 
estimate the size because of the effects of beam convolution and 
radial temperature gradients in the clumps. For this reason, and 
to check the general sensitivity to the LOS extent, we also cal¬ 
culated another set of models with FWHM at half of the value 
given above. Note that the clump size can also be significantly 
larger along LOS than in the plane of the sky. There is even 
some bias in this direction because the emission from elongated 
clumps and filaments becomes stronger when they are aligned 
along LOS. 

We performed the calculations iteratively, and on each itera¬ 
tion we solved the dust temperature distributions and calculating 
model predictions of the surface brightness that we then con¬ 
volved to the resolution of the observations. The calculations 
were carried out with a Monte Carlo radiative transfer program 
( Juvela|2005 1. The initial external radiation field corresponded to 
that of |Mathis et al.| ( p~983| >. The modeling was performed with 
two dust models. The first, in the following MWD, represents 
dust in normal diffuse interstellar medium ( |Li & Draine|2001~]l. 
The other one, in the following OH, was taken from Ossenkopf 
& Henning (1994 1 and corresponds to dust that has coagulated 
and accreted thin ice mantles over a period of 10 5 years at a den¬ 
sity 10 6 cm~ 3 . In our calculations, the dust opacities A:(250(rm) 
were 0.045 cm 2 g 1 for MWD and 0.22 cm 2 g 1 for OH, the 
value of Sect.[4]falling between the two. The value of k is crucial 
because it affects the column density, an important parameter of 
the subsequent line modeling. 

In the plane of the sky, the column density corresponding to 
each map pixel was adjusted by comparing the observed 350 pm 
surface brightness with the model prediction. The external radi¬ 
ation field was adjusted so that finally the 160 pm and 500 pm 
predictions also agreed with the observations to within 10 %. 
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Fig. A.2. Continued from Fig. 


A.l 


The observations could be fitted well with both dust models, 
with an external radiation field 70 - 80 % of the IMathis et al.l 
( |1983| values. Because we used background-subtracted surface 
brightness data, the model represents only the inner parts of the 
cloud without the diffuse envelope. Therefore, the radiation field 
in the models is somewhat weaker that the full radiation field 
outside the cloud. The column densities are listed in Table lC.il 
The values obtained with the OH dust model are close to the 
column densities estimated in Sect. [4] With MWD, the column 


densities are higher and in the northern point by more than the 
ratio of dust opacities, ~5. For a given dust opacity and radiation 
field, there is a maximum surface brightness that can be pro¬ 
duced with any column density. Thus, for a too small value of k, 
the column density of the model might be grossly overestimated. 
When estimated with the NICER method (using 2MASS stars, a 
spatial resolution of two arcminutes, and a value of Ry= 3.1), the 
visual extinction of northern clump is less than 3 mag. Taking 
the difference in the resolution into account, this is still consis- 
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tent with the Sect. [4] estimates and the result from the models 
with OH dust. However, the Ay measurement appears to rule 
out the values obtained with MWD dust (assuming the Ay is not 
severely underestimated either because of very clumpy column 
density structure or the presence of foreground stars) and give 
some support to the idea of dust opacity higher than that of dif¬ 
fuse medium. 


C.2. Modeling of the lJ CO(1-0) and C* 0(1-0) lines 

We modeled the 13 CO(1-0) and C 18 O(l-0) lines separately, tak¬ 
ing the density distribution directly from the continuum models. 
We only modeled the spectra toward the center of the clumps. 
The fractional abundance and kinetic temperature were assumed 
to be constant but in the non-LTE models the excitation temper¬ 
ature does vary and gives more weight to the dense regions. The 
velocity field was initialized by giving each cell a random veloc¬ 
ity vector with cr 3D = 1.0 km s 1 and assuming a turbulent line 
width with Doppler velocity Av D = V2cr 1D = 1.0 km s 1 . Thus 
the initial ’’microturbulence” within the cells is slightly larger 
than the ’’macroturbulence” between the cells. In the actual cal¬ 
culations, to match the observed line widths, both velocity com¬ 
ponents are scaled by the same number, typically by ~0.3 - 0.4, 
and the thermal line broadening is added to the turbulent line 
widths. The predictions of the models are optimized regarding 
the line width (the scaling mentioned above) and the line inten¬ 
sity that is adjusted by changing the molecular abundance. There 
are four cases for each position corresponding to two different 
assumptions of the LOS cloud extent and the two different dust 
models. Furthermore, some models were also calculated with 
a kinetic temperature of 7\j n = 12.0 K, instead of the default 
assumption of 7’ km = 10.0 K. Again, the true kinetic tempera¬ 
ture is unknown, but the comparison between 7’ k j n = 10.0 K and 
7 kin = 12.0 K gives some idea of the associated uncertainties. 

We performed the line calculations with the Monte Carlo 
program described in!|Juvela (1997 ». The obtained abundances 
are listed in Table Concentrating on the models based on 
continuum modeling with OH dust, the 13 CO values are on the 
order of the canonical value of 10 6 and the estimates are simi¬ 
lar for both assumed values of 7' km . However, because the abun¬ 
dances directly depend on the assumed column density, and thus 
the values of dust k, the absolute values are very uncertain. The 
relative abundance between l3 CO and C 18 0 should be more re¬ 
liable, as suggested by the similarity between the OH and MWD 
cases and the two cases of LOS density distribution. The abun¬ 
dance ratio is ~10 and thus larger than the value of 5.5 that was 
assumed in Sect. [3] 
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Table C.l. Result from radiative transfer modeling. The columns are (1) observed position, (2) column density estimated in Sect, [d] 
(3) width of the LOS density profile, (4) dust model, (5) column density of the continuum model, (6) assumed kinetic temperature, 
(7) 13 CO abundance of the model, and (8) relative abundance [ 13 CO]/[C ls O]. 


Position 

N(H 2 ) 

(10 21 cnr 2 ) 

LOS width 

Dust 

N(H 2 ) 

(10 21 cnr 2 ) 

7id n 

(K) 

[ lj CO] 
(10 6 ) 

[ 1J CO]/[C ls O] 

G131.65+9.75 N 

9.3 

Wide 

OH 

10.9 

10.0 

0.36 

11.4 

G131.65+9.75 N 

9.3 

Narrow 

OH 

12.3 

10.0 

0.32 

11.3 

G131.65+9.75 N 

9.3 

Wide 

MWD 

85.8 

10.0 

0.04 

10.4 

G131.65+9.75 N 

9.3 

Narrow 

MWD 

93.3 

10.0 

0.04 

10.5 

G131.65+9.75 S 

3.7 

Wide 

OH 

3.4 

10.0 

1.50 

16.2 

G131.65+9.75 S 

3.7 

Narrow 

OH 

3.7 

10.0 

1.22 

14.8 

G131.65+9.75 S 

3.7 

Wide 

MWD 

20.3 

10.0 

0.12 

10.7 

G131.65+9.75 S 

3.7 

Narrow 

MWD 

22.2 

10.0 

0.11 

10.8 

G131.65+9.75 N 

9.3 

Wide 

OH 

10.9 

12.0 

0.32 

10.7 

G131.65+9.75 N 

9.3 

Narrow 

OH 

12.3 

12.0 

0.29 

10.8 

G131.65+9.75 S 

3.7 

Wide 

OH 

3.4 

12.0 

1.11 

13.3 

G131.65+9.75 S 

3.7 

Narrow 

OH 

3.7 

12.0 

0.94 

12.7 
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Fig.B.l. The N 2 H + spectra observed toward G92.04+3.93, G93.21+9.55. G98.00+8.75, G105.57+10.39, and G132.12+8.95. The 
field G86.97-4.06 is shown in Fig. [4] The red dashed line is the fit to the hyperfine spectra and the lower subframe shows the fit 
residuals. The main beam temperature of the 23-12 component, the fitted radial velocity, the FWHM line width (km s 1 ), and the 
residual rms noise (as main beam temperature) are given in the frames. 
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